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O 
(N 

>>: 



in 



O 



(N 
> 
cn 
in 
o 
m 

o 



C. R. D' Angelo and a. L. Watts 

Instituut Anton Pannekoek, University of Amsterdam, Amsterdam 1098 XH, Tlie Netherlands 
Draft version March 2, 2013 

ABSTRACT 

Strong quasi-periodic oscillations (QPOs) in the tails of the giant gamma- ray flares seen in SGR 1806-20 and 
SGR 1900H-14 are thought to be produced by starquakes in the flaring magnetar. However, the large fractional 
amplitudes (up to ^^20%) observed are difficult to reconcile with predicted amplitudes of starquakes. Here 
we demonstrate that the steeply pulsed emission profile in the tail of the giant flare can enhance the observed 
amplitude of the underlying oscillation, analogously to a beam of light oscillating in and out of the line of sight. 
This mechanism will also broaden the feature in the power spectrum and introduce power at harmonics of the 
oscillation. The observed strength of the oscillation depends on the amplitude of the underlying starquake, the 
orientation and location of the emission on the surface of the star, and the gradient of the light curve profile. 
While the amplification of the signal can be significant, we demonstrate that, even with uncertainties in the 
emission geometry, this effect is not sufficient to produce the observed QPOs. This result excludes the direct 
observation of a starquake and suggests that the observed variations come from modulations in the intensity of 
the emission. 

Subject headings: pulsars: individual (SGR 1806-20, SGR 1900H-14), stars: magnetars, stars: oscillations, 
X-rays: stars 



1. INTRODUCTION 

Giant gamma-ray flares are thought to be produced by 
a large-scale rearra ngement of the super-strong ma gnetic 
field of a magnetar ( {Thompson & Duncanll 19951 120011) . The 
large amount of energy released during a flare (~10^^erg) 
likely triggers large-scale vibrations and quakes in the star's 
crust (Duncan 1998). The observation of strong, long- 
lived quasi-periodic oscillations (QPOs) in the decay tails 
of the giant flares of SGR 1900H-14 and SGR 1806-20 has 
been i nterpreted as evidence of such starquakes ( Israel et al] 
2005t IS trohmaver & Wattsl l2005t lmtts"& Strohmaver.,2006r 
Strohm aver & Watts 2006$ hereafter SW06), and much flieo- 



retical work has subsequently focused on the physics of star- 
quakes themselve s (e.g. .Levin 2006; Glampedakis et al. 200^ 
iGabler et al.ll201 U for a recent review see Watts 2011 ). 

Comparatively little work has focused on how the physical 
mo tion of the crust transl ates into a detectable QPO (although 
see lTimokhin et alll2008l) . This is not a trivial question: some 
of the QPOs have rms fractional amplitudes of up to 20% of 
the total flux in a pulse, whereas the maximum predicted am- 
plitude of a starquake (assuming the oscillation is restricted 
to the crust) is about O.Oli?* (where R^, is the stellar radius). 
A persistent starqu ake will more l ikely have an amplitude at 
least 10 X smaller (lDuncanlll998t iLevin & van Hovenll201 It 
M. Gabler, 2012 private communication). 

A solution could come from the steeply peaked pulse profile 
in the magnetar tail (see the pulse profile for SGR 1806-20 
in Figure j4|i. The observed QPOs are phase dependent, and 
are someti mes centered on the steep e d ge of the pulse (e.g.. 
Figure 3 of 'Strohmaver & Wattsll2005l) . Vstrohmaver & WattsI 
(|2005) suggested that beamed emission could enhance the ob- 
served strength of an underlying oscillation. As the beam 
edge sweeps through the observer's line of sight, the under- 
lying starquake will cause the edge to wiggle in and out of the 
observer's line of sight and amplify the signal. 

In this Letter we use a toy model for the emission to quan- 
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tify the effect that the pulse shape will have on an underly- 
ing small-amplitude oscillation. The observed oscillation de- 
pends strongly both on the pulse profile and its orientation on 
the surface of the star We demonstrate that the observed os- 
cillation can easily be several times larger than the underlying 
physical motion from a starquake. However, the enhancement 
is likely not enough to fully explain the observed amplitudes 
of the QPOs in SGR 1806-20 and SGR 1900H-14, suggest- 
ing that the observed oscillations are produced by a changing 
intensity in the emission rather than physical motion of the 
beam itself. 

2. MODEL FOR STARQUAKE 

To quantify how the signal from a starquake will be modi- 
fied by the pulse profile requires a description for the oscilla- 
tion and the emission. We adopt a simplistic model for both 
the emission and oscillation geometry, and further assume that 
the emission comes from close to the star's surface and moves 
with the crust. This can easily be generalized to a more de- 
tailed two-dimensional model for both the crust motion and 
the pulse pattern. 

Figure |T| shows a sketch of the star The star rotates with 
frequency fl^, about the z-axis, which is inclined by angle i 
from the observer's line of sight. The dashed line shows the 
observer's viewing parallel: rays parallel to the surface nor- 
mal on the viewing parallel will be observable. Assuming 
the emission is strongly beamed (necessary to produce sharp 
gradients), the emission from the star will be visible where 
it crosses the viewing parallel. We assume that the emission 
beam is a narrow strip in which the intensity varies strongly 
across the strip but stays roughly constant along it. The beam 
is seen in close-up in the figure and the line density shows the 
changing emission intensity. 

An oscillation with frequency i/q (where i^q ^ fi*) and 
physical displacement Ax is superposed on top of the star's 
overall rotation. The orientation can be random however we 
consider only the component perpendicular to the beam's gra- 
dient, since an oscillation along the gradient will not change 




Fig. 1. — Sketch of star with starquake. Main: the star rotates at f2*. with 
a starquake with amphtude Ax superposed on the rotation. The dashed line 
shows the observer's viewing parallel (at inclination i). The beam is observed 
(parallel hnes on the star's surface) where it crosses the viewing parallel. 
Inset: the viewing parallel (solid line) crosses the beam at angle a. A physical 
motion Ax is observed as a change in phase A(/>, which can be large if a is 
small. 

the emission profile. The beam crosses the viewing paral- 
lel with angle a, so the observed pulse profile is somewhat 
broadened. 

The high-frequency oscillation Ax corresponds to a shift in 
the rotation phase, A<i>: 



i?* sin i sin a 



sin(27ri'oO- 



(1) 



The 1 / sin a term arises from the misalignment between the 
beam and viewing parallel. As the beam oscillates in direction 
Ax, the change in flux corresponds to a shift in rotation phase 
A(j) larger than the physical displacement. This is because the 
observed beam is broadened along the viewing parallel, while 
Ax probes the 'true' beam profileQ. Angles a and i are un- 
known: however, we will argue later that they are insufficient 
to amplify a starquake into the observed features. 

Using the emission profile, P{^q), we can calculate 
and see how the starquake changes the overall variability in 
the light curve. Since A$ is small, observed profile is: 



P(<I>) = P($o + A$) 

= P($o) + A$f 



(2) 



OiA<P^ 



When the gradient of the light curve is large (~ 
P($o)/(A<I')), the starquake will be boosted and dominate 
the light curve. 

This process will have two other effects on the observed 
variability. First, since ^ changes, the observed amplitude of 
the oscillation varies, which will redistribute the signal power 
at into a range of frequencies. The width of the signal in 
the power spectrum can be estimated from the rate of change 

in the observed oscillation amplitude, i.e., This process 
will also introduce harmonics of i^o into the light curve. 
In summary: the strength of the oscillation seen by the 

viewer depends on Ax, i, a and | . | determines 
the width and shape of the feature in the power spectrum, as 
well as the magnitude of harmonics of i^o • 

' We are grateful to Yuii Levin for pointing this additional geometiical 
factor out to us. 



2.1. Making a Periodic Signal Quasi-Periodic 

The high-frequency signals detected in SGR 1806-20 and 
SGR 1900-1-14 had resolved widths suggesting that the under- 
lying oscillation is not a purely periodic signal. To understand 
how the pulse profile changes a signal with a range of frequen- 
cies, we add QPOs to the fight curve of SGR 1806-20. 

We model QPOs as noise processes with a Lorentzian dis- 
tribution of powers as input following the method suggested 
by .Timm er & Koenig (1995). The power spectrii m of the 
QPO in rms normalization (e.g., Ivan der KfisllT989h IS given 
by: 



PO 



<y^A /Ny 



2t: (A/2)2 + (j. - j.o)2 



(3) 



In Equation (O, /wy is the Nyquist frequency, ct^ is the QPO 
variance, A is its full-width at half-maximum (FWHM), and 
is its central frequency. 

The total power of the QPO is the same as for a sinusoid, 
only distributed over a range of frequencies. The variance of 
F{t) = ^0 sin(27ri'ot) over time interval T is: 



A2 fT a2 



(4) 



so that substituting the amplitude from Equation ([TJ, cr^ in 
Equation (O is given by: 



1 fAx 



1 



i?* sin i sin a 



(5) 



To produce a time series of length T, we generate a real- 
ization of its Fourier transform by taking the power (Equa- 
tion (|3]l) at each frequency, Vn = riAiy (Av = 1/T) and 
multiplying by a complex number {An + iBn) whose real and 
imaginary parts are drawn from a Gaussian random distribu- 
tion with a variance of one and mean of zero: 



fn — {Ar. 



iBn)\l ^5'QPo(i^n)- 



(6) 



The resulting time series has a variance and a mean of 
zero. This signal is the offset A$ to the star's phase $0 as a 
function of time. 

2.2. Example: A Power Law Light Curve 
We illustrate how the pulse profile changes an underly- 
ing oscillation with an emissionprofile P($o) « $q, adopt- 
ing different values of p. This idealized shape corresponds 
roughly to the edges of the pulse where QPOs have been ob- 
served. 

We use the same total flux for each profile, and neglect the 
effect of shot noise. We adopt a period of 10 s for the emission 
duration a time resolution of At = 9.5 x 10^^ s. We add 
a sinusoidal oscillation to the rotation phase, with amplitude 
Ax/R^, sin i sin a = 0.1 and frequency i^o = 600 Hz. 

Figure 12] shows the power spectra around 600 Hz, for f3 — 
[1,2,4,8], increasing from bottom to top . The inset of Figure 
|2] shows the different emission profiles. The rms variability 
from the starquake is measured by integrating the spectrum in 
an interval of 30 Hz centered at 600 Hz and subtracting the 
power contribution from the underlying emission profile. 

The main effect of steeper profiles is to increase the rms 
variabiHty of the starquake. For /3 = 1, the oscillation is a 
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Fig. 2. — Main: power spectra centered around iiigii-frequency feature for 
(bottom to top) /9 = [1, 2, 4, 8]. Tlie gradient in tlie liglit curve determines 
the width and amplitude of the resulting feature. Inset: the emission profile 
for each power spectrum. Solid:/3 = 1; dotted:/3 = 2; dashed:/? = 4; 
dot-dashed:^ = 8. 



TABLE 1 

Properties of Oscillation in Power-law Light Curve 
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I'D 


Width 


rms 


rms 




(Hz) 


(Hz) 


Fundamental 


1st harmonic 


1 


600 




2.2% 





2 


600 


0.11 


2.9% 


2.7 X 10"''% 


4 


600 


0.162 


8.4% 


1.2 X 10-3% 


8 


600 


0.3 


20.5% 


6.1 X 10-3% 



spike at 600 Hz with an rms amplitude of ^ 2%. As (3 in- 
creases, the ampHtude and width of the feature increase, so 
that for /3 = 8 the feature has an rms ampHtude of 20% — a 
ten-fold increase. 

For larger /?, the gradient of P($o) also changes substan- 
tially, which broadens the feature so that the FWHM for 
13 ~ 8 is 0.3 Hz. The full properties of the signal for each /3 
are summarized in Table [T| 

From the generated power spectrum we also measure the 
starquake harmonics. The last column of Table [T] shows the 
rms of the first harmonic (at 1200 Hz) for different /3. Even for 
/3 = 8, the power in the first harmonic is much smaller than 
the fundamental and would not be detectable above photon 
noise. 

These results show that a small oscillation can be substan- 
tially increased by a steep gradient of the light curve. We next 
apply the same technique to the light curve of SGR 1806-20, 
to demonstrate how the beam shape will affect the oscillation 
amplitude. 

3. THE TAIL OF SGR 1806-20 

The tail of the 2004 giant flai-e of SGR 1806-20 showed 
strongly pulsed emission with two main peaks in the 
pulse profile. QPOs appeared at different phases in the 
pulse profile, and also la sted for different lengths of time 
( IWatts & Strohmaveril2006i SW06). Two of the strongest 
QPOs, one at ~ 90 Hz and one at ~ 625 Hz are associated 
with steep gradients in the pulse profile. A dynamic power 
spectrum of the data showed that the 625 Hz QPO is strongest 



in the declining edge of the pulse (Figure 3 of SW06). We 
therefore focus our study on the same phase as in SW06. 

SW06 found phase-dependent QPOs at 93 Hz and 625 Hz 
in the SGR 1806-20 giant flare beginning 195 s after the main 
spike. The inset of Figure |4] shows the averaged pulse profile 
for the nine subsequent rotation cycles, and the solid lines 
show the phase segment when the QPOs were detected most 
strongly. 

A starquake will only be boosted by rotation-dependent 
variation of the light curve. The RXTE data of SGR 1806- 
20 shows substantial variability, particularly below 100 Hz, 
but the rotation period of the star is closer to 0.1 Hz. Ad- 
ditionally, at shorter time bins the signal-to-noise ratio drops 
and the variability from photon noise dominates. We thus bin 
and smooth the entire light curve to a time resolution of 0. 1 s, 
removing signal power at frequencies above ^ 10 Hz. 

Since the resulting feature depends on the gradient of the 
pulse profile and its derivative (Equation (|2]i), the binsize and 
smoothing change the detailed shape and intensity of the fea- 
ture. Experimenting with different binsizes and smoothing 
prescriptions shows the amplitude of the feature is robust (to 

^ 10 20%) to changes. The width and shape of the feature 

are more sensitive to higher derivatives so that the widths and 
the shapes of the features we produce serve more as qualita- 
tive examples than quantitative results. 

To probe the effect of the pulse shape on the ouput oscilla- 
tion strength, we make a light curve with two features, with 
A* = 0.1 sin(27ri^o t) at = 90 Hz and at 625 Hz. 
We then follow SW06 and examine the light curve 195 — 260 s 
after the main flare, calculating the power spectrum of the 3 s 
phase segment shown in the inset of Figure|4] The rms ampli- 
tude of the signal changes between rotations (from the chang- 
ing pulse profile) from 16% to 25%. Stacking and averaging 
the power spectra we measure the average power in each fea- 
ture. The stacked power spectrum centered around 90 Hz is 
shown in Figure |3] The feature has an amplitude of 21% and 
FWHM width of 1.6 Hz. The feature at 625 Hz has an rms 
amplitude of 20% and width 1.6 Hz. While the width is still 
smaller than the widths of the measured features (1.8 Hz for 
the 625 Hz feature and 8.8 Hz for the 93 Hz feature; SW06), 
it nonetheless demonstrates that a very narrow feature can be 
substantially broadened by the light curve profile. 

The harmonics of the features are stronger than for the 
power-law case with an rms amplitude of 3% and a large fre- 
quency spread in the first harmonic. This is still undetectable, 
even when the oscillation has an observed rms amplitude of 
20%. 

3.1. QPOs in SGR 1806-20 

The QPO features detected in SGR 1806-20 have measur- 
able widths between 1 17 Hz, and frequencies where pho- 
ton noise dominates the variability. Both these effects can 
change the measured properties of QPO, as we study below. 
We model the starquake as a quasi-periodic signal as out- 
lined in Section im defining an input width, central frequency 
and amplitude. To simulate photon noise we use the output 
light curve as a probability function to detect a photon in a 
given time bin and generate a light curve as a series of dis- 
crete events. Since these are both random noise processes, we 
generate a large number of realizations and measure the QPO 
properties statistically. 

As above, we add two QPOs to the phase of the SGR 1806- 
20 light curve and calculate the modified profile. One QPO is 
centered at 93 Hz with a width of 6 Hz, the second at 625 Hz 
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Fig. 3. — Feature produced by a sinusoidal oscillation with amplitude 0.1 
and frequency 90 Hz added to the phase for the light curve profile of SGR 
1806-20. 
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Fig. 4. — Main: power spectrum of a simulation using the SGR 1806-20 
light curve (smoothed to a resolution of 0.1s). QPOs were insetted at 93 Hz 
and 626 Hz. Inset: the rotation-averaged pulse profile 195-260 s after the 
main flare. The section of the phase profile used in the power spectrum is 
marked by vertical lines. 



with a width of 2 Hz. Both have an input amphtude of | A<i> | = 
0.1. We use the same rotation phase as above, stack the power 
spectra from nine successive rotations and bin the data with a 
binsize of 2.66 Hz. 

Figure |4] shows the power spectrum (in Leahy normaliza- 
tion) of one realization. The inset shows the averaged phase 
profile , with vertical lines denoting the phase segment used 
to make the power spectrum. This figure is comparable to 
Figure 2 of SW06. Both QPOs are visible by eye. In our 
spectrum the power above 10 Hz has been removed from 
the smoothing applied to the light curve. 

The QPO inserted at 93 Hz is detected with a mean width of 
6.6±1.5 Hz, rms amplitude of 21%±1.5%, and average cen- 
tral frequency of 9 1 .7±0.6 Hz. For the 626 Hz QPO, the mean 
width is 2.4±1.4 Hz, with an rms amphtude of 15%±1.5% 
and central frequency 624.9±0.35 Hz. 

The rms amplitude at 93 Hz is not significantly differ- 
ent from adding a sinusoid, but is measurably lower around 



626 Hz. In both cases, however, the strength of the feature 
is significantly stronger than the underlying oscillation. The 
features show modest broadening with a large dispersion be- 
tween different realizations. The systematic offset in the cen- 
tral frequency of the QPO appears because of the binsize. In 
simulations without photon noise, reducing the binsize shifts 
the central frequency of the QPO closer to the input one. Har- 
monics are not detected. 

4. DISCUSSION 

The steep pulse profile seen in the tail of giant magnetar 
flares substantially alters the observed properties of a star- 
quake, changing both its observed amplitude and power dis- 
tribution. While this effect is not large enough to transform a 
starquake with maximum amplitude A x/R^, ^0.001 — 0.01 
into a QPO with an rms amplitude of 20%, the boosting can 
still be considerable, especially if the emission comes from 
close to the rotation axis of the star, or if the steepest part 
of the beam gradient is misaligned with the viewing parallel 
(small a). 

There are several uncertainties in our calculation that could 
boost the amplitude of a starquake further The first of these is 
the unknown location of emission on the surface of the star. If 
the emission is close to the rotation axis, then a small physical 
motion produces a larger change in phase than at the equator 
However, to produce a ten-fold increase requires an inclina- 
tion angle of less than ^ 6°, which is improbable given the 
strong observed pulsations. 

The second way to boost the amplitude is for the beam to 
cross the observing parallel obliquely, as in Figure [T] Since 
in this case the actual beam edge is steeper than observed, 
the amplification is larger than expected from the observed 
beam shape. For the extreme case (a = 0) the beam runs 
exactly parallel to the line of sight so that no pulse is seen but 
there is a strong QPO. This scenario is unlikely, both that a 
is so small and that the beam edge is so much sharper (by a 
factor 1/1 sinaj) than observed. As an example: to produce 
a feature with a rms amplitude of 20% requires a gradient 
0.2 ~ dP/d<i> - ( Po) Aq (cf. Equation ©), which is - 
100 X steeper than the pulse profile in Figure |4] 

The input amplitudes required to produce the strong QPOs 
of SGR 1806-20 are large enough to produce QPOs that are 
observable throughout the entire profile. Applying the model 
across the entire pulse profile (using 3 s segments offset by 
O.IP* in phase), the oscillation has an output amplitude be- 
tween 10% and 25%-strong enough to be detected throughout 
the rotation. If the emission model is correct, this requires that 
the starquake amplitude varies significantly across the surface 
of the crust and the correlation with the falling pulse edge is 
somewhat coincidental. This is expected generically in star- 
quake models, however, and is seen directly in MHD simu- 
lations of coupled core-crust oscillations (Gabler et al. 201 1]) 
which show a strong phase dependence in the starquake am- 
plitude. 

The toy model suggests that the pulse profile can substan- 
tially influence the observed amplitude and other properties 
of a starquake. While we have used a simple model for the 
emission and starquake, this work can be easily generalized 
to use a physically-motivated model of the surface emission 
and starquake, in order to put constraints on the pattern of 
emission on the surface of the star during the quake. 

We thank Phil Uttley, Yuri Levin and Daniela Hup- 
penkothen for useful discussion and acknowledge support 
from an NWO Vidi Grant. 
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